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Light sterile neutrinos have been introduced as an explanation for a number of oscillation signals 
at Am 2 ~ 1 eV 2 . Neutrino oscillations at relatively short baselines provide a probe of these possible 
new states. This paper describes an accelerator-based experiment using neutral current coherent 
neutrino-nucleus scattering to strictly search for active-to-sterile neutrino oscillations. This exper- 
iment could, thus, definitively establish the existence of sterile neutrinos and provide constraints 
on their mixing parameters. A cyclotron-based proton beam can be directed to multiple targets, 
producing a low energy pion and muon decay-at-rest neutrino source with variable distance to a 
single detector. Two types of detectors are considered: a germanium-based detector inspired by the 
CDMS design and a liquid argon detector inspired by the proposed CLEAR experiment. 



I. INTRODUCTION 

Sterile neutrino models have been invoked to explain a 
series of intriguing oscillation signals at Am 2 ~ 1 eV 2 1- 
0]. These signals have relied on neutrino detection 
through charged current interactions. In the case of 
charged current appearance, the signal is interpreted as 
an active flavor oscillating to another active flavor, which 
can occur at these high Am 2 values if one or more neu- 
trino mass states with 777,4,... ~ I eV is added to the 
neutrino mass spectrum. The extra mass states are as- 
sumed to participate in neutrino oscillations, and must 
therefore be small admixtures of weakly-interacting neu- 
trino flavor states, with the remaining flavor composition 
being sterile (i.e. non- weakly-interacting) . In the case of 
charged current disappearance, the signal is interpreted 
as arising from active-flavor neutrino (e, fi, r) oscillation 
to any other neutrino flavor (e, /1, r, or s, with s being 
the sterile flavor). 

The oscillation probabilities for appearance and disap- 
pearance through charged current searches are expressed 
as functions of the active flavor content of the extra mass 
eigenstate(s) Q3I2]- In this paper, we assume that only 
one such extra mass state, 7714, exists. In that case, the 
oscillation probabilities are given by 



P(y a -> v 07ta ) = 4\U a4 \ 2 \U 04 \ 2 sin 2 {1.27 Am^L/E) 
in the case of active appearance searches, and 



(1) 



P{v a -> v*) = 4|C/ a4 | 2 (l - |f7 a4 | 2 )sin 2 (I.27Ar74L/£) 

in the case of active disappearance searches, where a, f3 = 
e, /it, t; fx corresponds to all flavors other than a, includ- 
ing active and sterile; |J7 Q 4| 2 corresponds to the a-flavor 
content of the fourth mass eigenstate; and L and E repre- 
sent the neutrino travel distance and energy, respectively. 
Note that neither search case is purely sensitive to the 
sterile neutrino content of the extra neutrino mass state, 
\U s4 \ 2 . In this paper, we discuss a strictly neutral current 



search using coherent neutrino scattering that allows for 
pure active-to-sterile oscillation sensitivity. 

Coherent neutrino-nucleus scattering is a well- 
predicted neutral current weak process with a high cross 
section in the standard model, as compared to other neu- 
trino interactions at similar energies. Despite this, the 
coherent interaction has never been observed as the keV- 
scale nuclear recoil signature is difficult to detect. The 
newest generation of ~10 keV threshold dark matter de- 
tectors provides sensitivity to coherent scattering [5] as 
the interaction signal is nearly identical to that which is 
expected from WIMP interactions. 

An active-to-sterile neutrino oscillation search is moti- 
vated in Section [TT] We describe an experimental design 
which makes use of a high intensity pion- and muon- 
decay-at-rest (DAR) neutrino source in Section III The 



coherent scattering process is introduced and the exper- 



imental design is discussed in Section IV Sensitivities 
to neutrino oscillations at Am 2 ~ I eV 2 are shown in 
Section El 



II. MOTIVATION FOR AN 
ACTIVE-TO-STERILE OSCILLATION SEARCH 

A decade ago, sterile neutrino oscillation models were 
largely motivated by the LSND anomaly \1\ 6 9 . This 
result presented a 3.8cr excess of v e events consistent with 
Vfx — > v e oscillations described by Eq. [I] at Am 2 ~ I eV 2 
and sin 2 20„ e = 4|[/ e4 | 2 |C/ M 4| 2 - 0.003. The apparent ap- 
pearance signal is thus interpreted as indirect evidence 
for at least one additional neutrino carrying the ability 
to mix with active flavors. Being mostly sterile, an ad- 
ditional neutrino avoids conflict with measurements of 
the Z invisible width [TUj (characteristic of three weakly- 
interacting light neutrino states) and the three-neutrino 
oscillation model established by solar [ITiri3] and atmo- 
spheric/accelerator [HHT7] experiments. 

The LSND signal was not present in a similar but less 



2 



sensitive — > D e oscillation search by the KARMEN ex- 
periment [18] ■ More recently, however, the MiniBooNE 
experiment [TH] has explored the Am 2 ~ 1 eV 2 parameter 
space and yielded a number of interesting results. Mini- 
BooNE features a higher beam energy and larger distance 
than LSND but preserves the L/E oscillation probabil- 
ity dependence, allowing for an independent cross check 
of the signal. In searching for v e appearance in a pure 
Vfj, beam, MiniBooNE has excluded — ► v e oscillations 
in the LSND Am 2 range at the 90% CL [20]. However, 
MiniBooNE's search for — > v e oscillations in "anti- 
neutrino-mode" is only consistent with the no oscillation 
hypothesis at the 0.5% level [2T]. The anti-neutrino re- 
sult is consistent with LSND and — > v e oscillations in 
the Am 2 = 0.1 — 1.0 eV 2 range. The statistics-limited 
measurement is expected to improve with additional data 
being taken through at least 2012. 

Recently, further results for v e disappearance at high 
Am 2 have been reported from short-baseline reactor anti- 
neutrino experiments. More specifically, a re-analysis of 
the anti-neutrino spectra produced by fission products in 
a reactor core [22] has led to an effect termed "the reactor 
anti-neutrino anomaly" , where the ratio of the observed 
anti-neutrino rate to the predicted rate deviates below 
unity at 98.6% CL [23 . This can be interpreted as dis- 
appearance according to Eq. [2] where charged current 
interactions of active flavors other than e are kinemati- 
cally forbidden, and/or where the oscillation was into a 
non-interacting sterile neutrino. Assuming CPT conser- 
vation, which requires that Eq.[2]holds for both neutrinos 
and anti-neutrinos, the strongest limits on D e disappear- 
ance come from a joint analysis of KARMEN and LSND 
v e + 12 C — > 12 N gs + e~ scattering events, analyzed for 
evidence of v e disappearance [24]. The reactor-anomaly 
signal is found to be marginally consistent with the KAR- 
MEN and LSND v e disappearance results. 

The above experiments feature a single source, sin- 
gle detector design. An alternative approach is a near- 
far detector configuration, where the measured flux in 
the near detector replaces the first-principles flux pre- 
diction. A near-far design removes a significant source 
of uncertainty due to the flux prediction, especially if 
the detectors are built to be nearly identical. Using the 
near-far technique, the CDHS [25], CCFR [25], and Sci- 
BooNE/MiniBooNE [27] experiments have probed neu- 
trino disappearance at Am 2 ~ 1 eV 2 using charged 
current interactions. Among the recent near-far compar- 
ison data sets, the MINOS experiment has set the only 
limits on active-to-sterile oscillations using neutral cur- 
rent interactions [2B]. The resulting limits using both 
charged current and neutral current interactions present 
a challenge in fitting sterile neutrino oscillation mod- 
els [29]. 

The aforementioned results underscore the experimen- 
tal and theoretical need for acquiring further data in ad- 
dressing the possibility of sterile neutrinos |30j . 



III. THE NEUTRINO SOURCE 

A DAR neutrino source can be employed to search for 
active-to-sterile neutrino oscillations through the neutral 
current coherent scattering interaction. DAR neutrinos 
have been identified as an excellent source for neutrino- 
nucleus coherent scattering studies [3"TH3"3"] because the 
neutrinos are produced in an energy region (<52.8 MeV) 
where the coherent neutrino scattering cross section is 
higher than all others by about one order of magnitude. 

A search for active-to-sterile oscillations is envisioned 
with a series of measurements at different values of L 
from the DAR source. In our design, a cyclotron directs 
a proton beam to two graphite targets embedded in a 
single iron shield. As the DAR neutrino flavor content 
and energy distribution is driven by the weak interaction, 
the well understood flux emitted isotropically from each 
target will be effectively identical, barring oscillations, at 
each baseline L. 

As discussed in Ref. [5J, an ideal neutrino interaction 
target for a DAR source is a direct dark matter detection 
device sensitive to keV-scale nuclear recoils. We consider 
two dark matter detector technologies; a germanium- 
based CDMS-style detector [M] and a liquid argon-based 
one similar to the CLEAN [35] and CLEAR [33] designs. 

A. About decay-at-rest sources 

In the DAR source described here, neutrino production 
begins with 800 MeV protons impinging on a target to 
produce low energy charged pions primarily through the 
A resonance decay. The pion decay chain tt + — > n + Vu, 
/i + — » e + D^v e , produces the neutrino flux shown in Fig.nl 
The v e content that arises from iv~ production and sut> 
sequent decay-in-flight (DIF) is well below 10 -3 [5j [TS] 
due to 7T _ capture. These features provide an ideal 
source for neutrino appearance searches [6j [TH] EU [37] , 
active-to-sterile searches relying on the charged current 
interaction [23], and an active-to-sterile neutral current 
search [35] as discussed in this paper. 

B. Targeting to allow multiple baselines 

A high intensity source of 800 MeV protons is being 
developed by the DAE5ALUS collaboration [37]. This 
design utilizes cyclotron-based accelerators [35] 00] in- 
stalled at three sites near a very large water- or oil-based 
detector. The experiment described here could use one 
of these DAEDALUS cyclotrons combined with a dual- 
target configuration as a neutrino source. 

In a baseline scenario, the cyclotron-based beam will 
be diverted between the two targets so that no target re- 
ceives more than 1 MW average power. The beam will 
be directed at 90° with respect to the detector, so as to 
minimize DIF backgrounds. Notably, a multi-target de- 
sign can also be exploited for a charged current neutrino 
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interaction oscillation measurement with a common de- 
tector and multiple baselines. The main technical issue 
in the two-target cyclotron design is maintaining a good 
vacuum in the two-prong extraction line. The beam will 
be "painted" across the face of each target in order to 
prevent hot spots in the graphite, an effect which will 
dominate the ±25 cm uncertainty on the experimental L 
from each neutrino source. The targets will be arranged 
in a row enveloped within a single iron shield, with the 
detector located 20 m downstream of the near target and 
40 m downstream of the far target. This configuration 
has been found to provide the best overall sensitivity to 
the LSND allowed region. 

The analysis below exploits the L dependence of neu- 
trino oscillations. Therefore, the flux of protons on each 
target must be well understood in time; standard proton 
beam monitors allow a 0.5% measurement precision. The 
absolute neutrino flux is less important, as sensitivity to 
the oscillation signal depends on relative detected rates 
at the various distances. The systematic uncertainty as- 
sociated with the flux normalization is 10% if there is no 
large water or oil detector available and 1.1% if such a 
detector does exist [36 . A high statistics i/-electron scat- 
tering measurement at a large water detector provides a 
precise determination of the flux normalization. 

IV. DETECTING COHERENT NEUTRINO 
SCATTERING 



the number of neutrons, number of protons, and weak 
mixing angle, respectively]; M is the nuclear target mass; 
T is the nuclear recoil energy; and E v is the incoming 
neutrino energy. The ~5% cross section uncertainty, the 
actual value depending on the particular nuclear target 
employed, is dominated by the form factor |42j . 

Coherent neutrino scattering is relevant for the under- 
standing of type II supernova evolution and the future de- 
scription of terrestrial supernova neutrino spectra. Mea- 
suring the cross section of the process also provides sensi- 
tivity to non-standard neutrino interactions (NSI) and a 
sin 2 Qw measurement at low Q [31] . Cross section mea- 
surements as a function of energy on multiple nuclear 
targets can allow the cross section dependence on NSI 
and 9\v to be isolated and understood. As demonstrated 
here, neutrino oscillations can also be cleanly probed us- 
ing coherent scattering. 

The difficulty of coherent neutrino scattering detection 
arises from the extremely low energy of the nuclear recoil 
signature. For example, a 20 MeV neutrino produces a 
maximum recoil energy of about 21 keV when scattering 
on argon. Both a CDMS-style germanium detector [33] 
and a single phase liquid argon detector, such as the one 
proposed for the CLEAR experiment [33], are consid- 
ered in this paper for detecting these low energy events. 
Other dark matter style detector technologies, especially 
those with ultra-low energy thresholds, can be effective 
for studying coherent neutrino scattering as well. 



Coherent neutrino-nucleus scattering, in which an in- 
coming neutrino scatters off an entire nucleus via neu- 
tral current Z exchange [41] . has never been observed 
despite its well predicted and comparatively large stan- 
dard model cross section. The coherent scattering cross 
section is 
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where Gp is the Fermi constant; Qw is the weak charge 
[Q w = N - (1 - 4 sin 2 ^)^, with N, Z, and 9 W as 
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FIG. 1: Energy distribution of neutrinos from a DAR source. 



A. Experimental Setup 

The envisioned experimental setup is consistent with 
the current DAEDALUS accelerator proposal and follows 
a realistic detector design. A single DAEDALUS cy- 
clotron will produce 4 x 10 22 i//flavor/year running with 
a duty cycle between 13% and 20% (33139]. A duty cy- 
cle of 13% and a physics run exposure of five total years 
are assumed here. With baselines of 20 m and 40 m, 
the beam time exposure distribution at the two baselines 
is optimal in a 1 : 4 ratio: one cycle to near (20 m), 
four cycles to far (40 m). Instantaneous cycling between 
targets is important for target cooling and removes sys- 
tematics between near and far baselines associated with 
detector changes over time. The accelerator and detector 
location is envisioned inside an adit leading into a sharp 
300 ft rise at the Sanford Research Facility at Homes- 
take, in South Dakota. The neutrino flux normalization 
uncertainty at each baseline is conservatively expected 
at 1.5%. We assume the flux has been constrained to 
this level by an independent measurement of i/-electron 
scattering with a large water-based Cerenkov detector 
also assumed to be in operation at Sanford Labs. The 
1.5% uncertainty estimate takes into consideration the 
theoretical uncertainty in the ^-electron scattering cross 
section and the statistics achievable with a large water 
detector. The flux normalization correlation coefficient 
between the near and far baselines is conservatively set 
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v source 


4 X 10^ i//fiavor/ycar 


1-JUKiy IdLliUl 


±o/o 


Baseline correlation 


0.99 


v flux norm, uncertainty 


1.5% 


Uncorr. sys. uncertainty 


0.5% 


Distances from v source 


20 m, 40 m 


Exposure 


5 years: 1 near, 4 far 


Depth 


300 ft 



TABLE I: The experimental configuration assumptions. 



to 0.99, its deviation from unity being dominated by dif- 
ferences between the two beam dumps. An uncorrelated 
systematic uncertainty of 0.5% at each baseline, is also 
included. The general experimental assumptions can be 
seen in Table |U 

We also consider a "dedicated" physics run scenario in 
which the duty factor is raised from 13% to 50% for all 
five years. With the instantaneous power achievable re- 
maining constant, this change leads to an average power 
increase of a factor of 4. Steady-state and beam-related 
backgrounds also increase by this factor in a dedicated 
scenario. The nominal duty factor of 13% is driven by 
the requirement that the various DAEdALUS accelerator 
baseline beam windows do not overlap in time. A dedi- 
cated scenario is possible in consideration of maintaining 
sufficient target cooling and the phased DAEDALUS de- 
ployment timeline. The timeline calls for a cyclotron or 
set of cyclotrons installed exclusively at a single "near" 
baseline, close to a large water detector, for at least 
five years [37] ■ With a 13% duty factor only required 
when all baselines have operational accelerators, a longer 
duty factor and higher average power seems possible in 
DAE<5ALUS single-baseline-only operation. Note that al- 
though only two targets are required for the experimen- 
tal design described here, supplementing the bcamlinc 
with more targets can ensure optimal use of beam time 
in consideration of cooling requirements and ultimately 
increase neutrino oscillation sensitivity. 



1. Germanium detector - signal and backgrounds 

A low-threshold germanium-based detector, such as 
CDMS, measures phonons and ionization from electronic 
and nuclear recoils 43J. A CDMS detector consists of a 
large germanium crystal (0.25 — 1 kg) operated at cryo- 
genic temperatures (^100 mK) with thousands of super- 
conducting transition-edge sensors (TESs) photolitho- 
graphically patterned on the top and bottom surfaces. 
The TESs are wired in parallel to form four readout chan- 
nels on each surface, which measure phonons created in 
particle interactions. The particle-induced ionization is 
also measured by electrodes on the crystal surface. The 
ratio of the energy in these two channels is a powerful 
discriminator between nuclear and electronic recoils. 

A 100 kg active mass of germanium is considered for 
the experiment described here, similar to proposed dark 



matter searches [23]. The detection efficiency above a 
10 keV threshold is set to 0.67 with a 3% energy reso- 
lution near the threshold. These assumptions are rea- 
sonably conservative and consistent with future expecta- 
tions [23 131 • 

Two classes of background events are considered for a 
germanium detector: 

1. Misidentified electronic recoils - Electronic recoils 
can be produced by photons and beta parti- 
cles interacting with the active detection medium. 
Misidentification of such events is particularly 
problematic near the detector surfaces, where 
the collection of electron-hole pairs is suppressed 
and discrimination is less effective. Existing ex- 
periments have demonstrated an electronic recoil 
misidentification rate of less than 1 event per 
100 kg-days exposure [34]. Upgrades to detector 
design are expected to improve discrimination by a 
factor of 10 4 [2S]. The assumed rate of radiogenic 
background detection (~2 events/year) is negligi- 
ble. 

2. Cosmogenic neutrons - Single scatter neutrons can 
produce a signal identical to a coherent neutrino 
scattering event, and the rate of these events would 
be significant at a shallow site. As a point of ref- 
erence for surface experiments, the CDMS exper- 
iment located at the Stanford Underground Fa- 
cility with 16 m.w.e. of overburden measured a 
neutron background of 0.67 events/ (kg- day) [47] . 
This figure could be significantly reduced with 
additional active and passive shielding and the 
larger overburden envisioned for the DAE5ALUS 
site. A cosmogenic-induced background of 0.1 de- 
tected events/(10 kg-day), after correcting for effi- 
ciency and during beam-on, is assumed. This value 
is considered a design goal and can be met with a 
300 ft overburden and modest active and passive 
shielding. 

In this study, the estimated radiogenic and cosmogenic 
background rates are distributed evenly across the ger- 
manium nuclear recoil energy range considered, 10 keV 
to 100 keV. 



2. Liquid argon detector - signal and backgrounds 

A single phase liquid argon detector can be used to de- 
tect the scintillation light created by WIMP- or coherent 
neutrino-induced nuclear recoils. Such detectors employ 
a large, homogeneous liquid argon volume surrounded by 
photomultiplier tubes (PMTs). Inner detector surfaces 
as well as the PMTs themselves are usually covered in a 
wavelength shifting substance which converts the 128 nm 
scintillation light into the visible spectrum for detection. 

A 456 kg active mass of liquid argon with a flat effi- 
ciency of 0.50 above a 30 keV energy threshold is con- 
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sidered for the experiment described here. The detec- 
tion volume and efficiency are consistent with the pro- 
posed CLEAR design [33]. An 18% energy resolution 
near threshold is used, assuming resolution slightly worse 
than what would be expected from photoelectron Pois- 
son statistics [48], 6 photoelectrons/keVee light collec- 
tion, and a quenching factor of 0.25 |49j . 

There are three primary sources of background that 
are considered for a single phase liquid argon detector: 

1. Cosmogenic neutrons - The muon-induced neu- 
tron background is contingent on the geometry of 
the site, overburden, and active/passive detector 
shielding. Muon events and muon-induced neu- 
trons can be vetoed with high efficiency and low de- 
tector dead time in a liquid argon detector near the 
surface [33]. The target design cosmogenic back- 
ground is 0.1 detected events/(10 kg-day). 

2. 39 Ar contamination - 39 Ar is a naturally occurring 
radioisotope with an isotopic abundance of 39 Ar/Ar 
= 8 x 10~ 16 , corresponding to a specific activity of 
1.01 Bq/kg [SD]. The isotope is a beta emitter with 
an energy endpoint of 565 keV. 

Pulse-shape discrimination (PSD) can be used to 
separate the 39 Ar-induced electronic recoils from 
the nuclear recoils produced in WIMP and coher- 
ent neutrino scattering events |51j . The electronic 
recoil contamination (ERC) of the nuclear recoils 
decreases exponentially as the number of photo- 
electrons detected increases. Reference [52] mea- 
sures the ERC for PSD in the single phase liq- 
uid argon DEAP-1 detector (4.85 photoelectrons 
per keVee) and also provides a "theoretical" Monte 
Carlo estimate of the ERC attainable for an ideal 
detector with 4-7T PMT coverage and 6 photoelec- 
trons/keVee. Both scenarios correspond to a 50% 
efficiency for nuclear recoil detection in the fiducial 
volume. Note that, according to Ref. [S3], the Mi- 
croCLEAN experiment has achieved 6 photoelec- 
trons/keVee sensitivity. The abundant 39 Ar back- 
ground could further be alleviated with the use of 
depleted argon from underground sources, which 
has an isotopic abundance of 39 Ar that is <5% of 
natural argon at the surface [54] . Figure [2] shows 
the rate of 39 Ar after PSD with 13% on-time, for 
two assumptions of ERC reported in Ref. [52] . The 
theoretical ERC with non-depleted liquid argon is 
employed for this study. 

3. Surface contamination - Radioactive impurities on 
the detector surfaces can decay and contribute to 
the background. These surface backgrounds have 
been measured in the DEAP-1 detector and were 
found to have an activity of 1.3 x 10~ 4 Bq [55J. De- 
pending on the origin of these events, the scaling 
and resulting background prediction can differ. If 
the events are due exclusively to 210 Pb surface con- 
tamination, the DEAP-1 figure can be scaled by de- 




Measured (4.85 pe/keVee) 
Theoretical (6 pe/keVee) 



Recoil energy (keVr) 

FIG. 2: The expected 39 Ar background energy spectrum un- 
der two sets of assumptions. The line labeled "Measured" 
corresponds to an ERC that was obtained in a detector with 
4.85 photoelectrons/keVee. The line labeled "Theoretical" 
is the ERC simulated in an ideal detector with 6 photoelec- 
trons/keVee and represents the background ERC used for this 
study. Both lines correspond to a 50% efficiency for detecting 
nuclear recoils in the fiducial volume. 



tector surface area to yield 1.3 x 10 4 events/year. 
However, this rate may be substantially reduced 
by the use of cleaner materials, scrubbing of the 
surface, and fiducialization. A surface background 
contamination of 100 detected events/year is as- 
sumed here and can be considered a design goal. 

The 30 keV energy threshold employed here is larger 
than the oft-chosen 20 keV threshold in single phase liq- 
uid argon detectors in order to mitigate the steeply falling 
39 Ar contamination. If 39 Ar discrimination improves in 
a future design, adjusting the threshold to (e.g.) 20 keV 
would allow a 60% larger signal sample. In this study, the 
estimated surface and cosmogenic background rates are 
distributed evenly across the argon nucleus recoil energy 
range considered, 30 keV to 200 keV. 

One additional possibility that would significantly re- 
duce the non-beam-related background would be to use a 
pulsed source of neutrinos, such as at the Spallation Neu- 
tron Source (SNS). The SNS produces protons in very 
short bunches of <750 nsec at a rate of about 60 Hz, 
so that the time window for expected signal events is a 
small fraction of the total running time. Combining a 
pulsed DAR beam structure with a liquid argon detector 
was previously proposed by the CLEAR experiment [33j , 
allowing them to claim an additional rejection of 6 x 10 -4 
for steady-state, non-beam-related backgrounds using a 
timing cut. Although the DAE<5ALUS proposal does not 
include this timing structure, the experimental concept 
described here could be employed at other facilities. 

The detector-specific assumptions are summarized in 



() 





'°Ge 


4U Ar 


Active mass 
Efficiency 
Threshold 
at threshold 
Radiogenic background 
Cosmogenic background 
Beam-related background 


100 kg 
0.67 (flat) 
10 keV 

3% 
2 / year 
0.1/(10 kg-day) 
/ year 


456 kg 
0.50 (flat) 
30 keV 

18% 
See text 
0.1/(10 kg-day) 
0/year 



TABLE II: The assumptions relevant for the specific detector 
technologies considered. 



Table |TT] and the expected signal and background rates 
are shown in Fig. [3j 




20 40 60 80 100 120 140 160 180 200 

Recoil energy (keVr) 

FIG. 3: The expected non-oscillated signal and total back- 
ground rates at a 20 m baseline for the two detector tech- 
nologies considered in the baseline physics run scenario. The 
rates correspond to what is expected for one full year of near- 
target-only operation at 13% duty factor. 



tector far enough from the source that the beam-related 
neutron flux is negligible. 

A precise estimate of the neutron flux would require 
detailed knowledge of the experimental site, beam con- 
figuration, and shielding. The neutron flux is estimated 
with a Monte Carlo simulation of the experimental geom- 
etry consistent with the DAEDALUS proposal [37] and 
several simplifying assumptions. Instead of simulating 
the passage of neutrons through the beam dump shield- 
ing, we simply assume that a cubic shield with sides of 
length 6 m is sufficient to reduce the escaping neutron 
flux to a level consistent with safety regulations. Also, 
we assume that this cube of shielding is adjacent to a 
rock (Si02) cliff. The maximum permissible annual dose 
for workers in a restricted area with a neutron beam is 
100 mRem 56 . The neutron flux escaping the shielding 
is set to a rate equivalent to an exposure of 100 mRem 
in 40 hours. 

Using a Geant4-based simulation [57J , neutrons are in- 
jected at the edge of the shielding cube. The neutrons 
are simulated in energy bins from 0-30 MeV. The flux 
is tallied at 20 cm intervals into the rock cliff, and the 
fluxes beyond 1 m into the cliff are fit to the functional 
form 



where A and A are fit parameters, and z is the distance 
from each flux tally point to the DAR source. The neu- 
tron fluxes are in reasonable agreement with this func- 
tional form. The fit function is then used to extrapolate 
the flux to a full year of running and larger distances from 
the source. A simulation is also employed to estimate the 
fraction of incident neutrons that produce single-scatter 
nuclear recoils in the detection volume. Less than 0.2 
beam-related events are expected per year for a 456 kg 
liquid argon detector at a 12 m baseline. The beam- 
related background at 20 m from the source, the shortest 
relevant detector baseline considered here, is therefore 
assumed to be negligible. 



3. Neutron flux from the source 

DAR sources produce a large flux of neutrons, aris- 
ing from spallation reactions of protons with the beam 
dump material. For the DAR source considered here, 
the neutrons have energies up to 800 MeV. In a 1 MW 
beam, the neutron production rate is ^10 16 /s and at- 
tenuation lengths may be as high as tens of centimeters. 
Single scatter neutrons can produce elastic recoils in the 
detector volume that are indistinguishable from coherent 
neutrino scattering on an event-by-event basis. More- 
over, because the neutron flux is attenuated by matter, 
underestimating the neutron background in the detec- 
tor could mimic a deviation from the l/r 2 -dependence of 
the neutrino flux, similar to what is expected for neutrino 
disappearance. It is therefore essential to locate the de- 



V. MEASUREMENT STRATEGY AND 
SENSITIVITY 

A. Overall strategy 

Neutrino oscillations depend upon neutrino energy and 
distance traveled. Since the neutrino energy cannot be 
reconstructed precisely with the coherent interaction, our 
sensitivity to the oscillatory behavior arises mainly from 
L, a value which is well determined by the location of the 
target being used at any given time and its distance to 
the common detector. In the case that a disappearance 
signal is detected, the target exposure priorities for the 
two baselines can be optimized to maximize sensitivity. 

The purely neutral current experiment described is 
sensitive to the effective disappearance of all three types 
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of neutrinos present in the beam, v^, z? M , and v e , into 
v s . We assume this disappearance can be approximated 
by a two-neutrino oscillation driven by a Am 2 in the 
LSND allowed region, and that the oscillation proba- 
bility under the approximation is the same for neutri- 
nos and anti-neutrinos. The baselines for the experi- 
ment, 20 m and 40 m, have been chosen in order to 
provide the best sensitivity to the LSND allowed param- 
eter space, given the neutrino energy spectrum of each 
flavor in the beam. The experiment described here pro- 
vides indirect sensitivity to the LSND allowed param- 
eter space by simultaneously measuring terms describ- 
ing the amplitude of active neutrino mixing to a sterile 
neutrino: 4|[/ e 4| 2 |t/ S 4| 2 in the case of u e in the beam, 
and 4|£/ /J 4| 2 |[/ S 4| 2 in the case of za, and in the beam. 
These terms are then translated to the appearance am- 
plitude measured by LSND, sin 2 20 Me = 4|J7 e4 | 2 |lV 4 | 2 - 
Sensitivity to sin 2 29 ^ e , along with simultaneous sensi- 
tivity to sin 2 26» ee = 4|L/ e4 | 2 (l - \U ei \ 2 ) and sin 2 29^ = 
4| 1 2 (1 — 1 1 2 ) , are considered the figures of merit 
here, as they can be easily compared to existing charged 
current appearance and disappearance measurements. 
Of course, distinguishing between sin 2 26 ee and sin 2 29 ^ 
in the case of an observed disappearance is not possi- 
ble in a flavor-blind experiment. Therefore, we rely on 
marginalizing over the full parameter space of {U^l and 
\U e 4\ explored, in the most conservative case possible, 
when drawing sensitivity contours for each case. 

The sensitivity to any particular set of oscillation pa- 
rameters is obtained by simultaneously fitting the ex- 
pected flavor-summed coherent signal events as a func- 
tion of recoil energy at the near and far baselines. The 
events at each baseline are distributed among bins of nu- 
clear recoil energy (1 bin/10 keV); however, the sensitiv- 
ity results are largely insensitive to the number of recoil 
energy bins used in the comparison. 

B. Sensitivities 

The signal predictions are evaluated for each set of os- 
cillation parameters, Amjj = Am 2 , | C/ M 4 1 , and \U e 4\. A 
X 2 is calculated by comparing the oscillations-predicted 
spectra, including backgrounds, to the no-oscillations 
prediction. 

The x 2 is constructed as 

where i and j denote the energy bins at the near and 
far baselines, respectively; Pi is the oscillations-predicted 
event spectrum as a function of N h [ ns = 1, 10, 11, 20 
bins, corresponding to (e.g.) 10 energy bins for the two 
baselines appended side by side; N is the correspond- 
ing no-oscillations spectrum; and M~^ 1 is the inverse co- 
variance matrix including statistical and systematic un- 
certainties and normalization systematic correlations be- 



tween the two baselines and different recoil energy bins. 
Note that the background contributions to Pi and Ni can- 
cel. The background-contributed statistical uncertainty, 
however, is accounted for in M^ . The background contri- 
bution can be measured with high statistics during beam- 
off cycles, and so systematic uncertainties associated with 
background are small relative to statistical uncertainties. 

The oscillations-predicted spectra, Pi , are obtained by 
summing over all neutrino flavors predicted in each recoil 
energy bin of the unoscillated spectrum, and reweighting 
each neutrino according to its flavor a — e,/i by the 
following "active" survival probability 

P(v a -> Active) = 1 - P{v a V s ) 

= 1 - sin 2 20 as sin 2 (1.27Am 2 L/£;)(6) 

where inactive can be any active state including v ai and 
sin 2 29 as — 4|£/ Q 4| 2 |£/ S 4| 2 . By unitarity assumptions, 
|[/ s4 | 2 is a function of J2 a =e^,r \ u a4\ 2 , 

\U si \ 2 = l- £ \U a4 \ 2 . (7) 

a=e,fj,,T 

During the fit, we vary |/7 e 4|, {U^], and Am 2 . For sim- 
plicity, however, we assume \U T 4\ = 0. Note that a non- 
zero {Ur^l would increase the active survival probability 
for any given \U e i\ and |L^i4|, and would therefore make 
this search slightly less sensitive to oscillations in terms 
of sin 2 26> Me . On the other hand, if non-zero \U e 4\ and 
| U^4 1 were to be established independently by other short 
baseline experiments, the type of neutral current search 
outlined in this paper may offer sensitivity to U T 4, de- 
pending on the sizes of \U e 4,\, {U^l and \U T 4\. 

Figures |2J [5] and [6] show the expected sensitivity to 
the LSND allowed region with a germanium detector in 
the baseline and dedicated physics run scenarios and an 
argon detector in the baseline scenario, respectively. In 
obtaining the sensitivity curves, the 3D search grid is 
reduced from (Am 2 , |?7 e 4| 2 , |£^ t 4| 2 ) to a 2D space of 
Am 2 and sin 2 29^ = 4|t/ e4 | 2 |t/ At 4| 2 . Note that a non- 
zero sin 2 29fj, e requires both v e and disappearance. 

The sin 2 29^ e sensitivity curves are obtained using a 
raster scan in Am 2 space. That is, each curve maps out 
the maximum sin 2 29^ e = 4|[/ e4 | 2 |l/ Al4 | 2 which satisfies 
X 2 < Ax 2 u t at a given confidence level, for each point 
in Am 2 . The 90%, 99%, and 3cr confidence level curves 
shown in this paper correspond to Ax\ ut =1.64, 6.63, 
and 9.00 for a one degree of freedom, one-sided raster 
scan (90%), and a one degree of freedom, two-sided raster 
scan (99% and 3<r), respectively. 

Figure [7] shows the oscillation sensitivity for a germa- 
nium detector in terms of the disappearance amplitudes 
which would be accessible in charged current searches, 
sin 2 29 ee = 4|[/ e4 | 2 (l-|[/ e 4| 2 ) and sin 2 29^ = 4|[/„ 4 | 2 (1- 
| 1 2 ) overlaid with the region allowed by LSND at 
90% CL, assuming the LSND best-fit Am 2 = 1.2 eV 2 . 
The curves are obtained using a one-sided raster scan 
in sin 2 29 ee with the Ax 2 ut values defined above. The 
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FIG. 4: Sensitivity to the LSND 90% CL allowed parameter 
space with a germanium-based detector under the baseline 
physics run scenario. 
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FIG. 6: Sensitivity to the LSND 90% CL allowed parameter 
space with an argon-based detector under the baseline physics 
run scenario. 
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FIG. 5: Sensitivity to the LSND 90% CL allowed parameter 
space with a germanium-based detector under the dedicated 
physics run scenario. 



current oscillation search correspond to sin 2 29^ < 0.1 
at 90% CL, for Am 2 = 1.2 eV 2 [28]. 

Figures [3] and [6] show that, despite the difference in 
fiducial mass, the 100 kg germanium detector performs 
slightly better than the 456 kg liquid argon one. The 
difference is in part due to the difference in nuclear re- 
coil energy threshold; 10 keV for germanium, 30 keV 
for argon. This emphasizes the fact that a low detec- 
tor energy threshold is important for obtaining a high- 
statistics sample of coherent neutrino scattering events 
as the rate is dominated by events with very low energy 
recoils (<10 keV). 

In a baseline physics run scenario, an experiment fea- 
turing a germanium- or argon-based detector can exclude 
the LSND best-fit mass splitting (Am 2 = 1.2 eV 2 ) at 
3.8ct or 3.4ct, respectively. The LSND best-fit mass split- 
ting is excluded at 4.8cr in the dedicated, germanium- 
based physics run scenario considered. For sensitivity in 
terms of sin 2 29 ee and sin 2 29^ , a germanium-based ex- 
periment in the baseline scenario could exclude nearly all 
of the available 90% CL LSND parameter space at the 3<7 
level and large portions of the available reactor anomaly 
allowed region, assuming Am 2 ~ 1.2 eV 2 . 



figure also shows the approximate region of sin 2 29 ee val- 
ues allowed at 90% CL by fits to the reactor anomaly 
and gallium experiment calibration data sets in Ref. [23] . 
The "reactor" allowed contour is for Am 2 >1.5 eV 2 and 
is relatively independent of Am 2 in this region. As a 
reference, limits on sin 2 26* MM from the MINOS neutral- 



VI. CONCLUSIONS 

This paper has described a method to search for 
active-to-sterile neutrino oscillations at relatively short 
baselines using neutral current coherent neutrino-nucleus 
scattering. Detection of such a process could definitively 
establish the existence of sterile neutrinos and measure 
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FIG. 7: Sensitivity to disappearance amplitudes accessible 
with charged current searches, assuming the LSND best-fit 
Am 2 — 1.2 eV 2 . The sensitivity corresponds to a germanium- 
based detector under the baseline physics run scenario. The 
LSND band represents the 90% CL allowed values of sin 2 28^ e 
at Am 2 = 1.2 eV 2 . "Reactor" refers to the result reported in 
Ref. [23] and indicates the range of sin 2 26 ee values preferred 
by a joint fit to reactor and gallium experiment calibration 
measurements. The reactor result is nearly independent of 
Am 2 , for Am 2 values near and above 1.5 eV 2 . 



their mixing parameters. 

An experiment that relies on the high statistics de- 
tection of an as-yet-undetected process is obviously dif- 
ficult. However, all of the technology required for such 
an experiment either exists or has been proposed with 
realistic assumptions. A cyclotron-based proton beam 
can be directed to a set of targets, producing a low en- 
ergy neutrino source with multiple baselines. This allows 
a measurement of the distance dependence of an oscil- 
lation signal without moving detectors or instrumenting 
multiple devices. Both a germanium-based detector in- 
spired by the CDMS design and a liquid argon detector 
inspired by the proposed CLEAR experiment would be 
effective for performing these measurements. 

Along with relevance in understanding Type II super- 
nova evolution and supernova neutrino detection, coher- 
ent neutrino-nucleus scattering can provide sensitivity to 
non-standard neutrino interactions, the weak mixing an- 
gle, and, as shown in this paper, neutrino oscillations 
at Am 2 ~ 1 eV 2 . Depending on the detector technology 
and run scenario, the experiment described is sensitive to 
the LSND best-fit mass splitting at the level of 3-5cr and 
can probe large regions of the LSND and reactor anomaly 
allowed regions. The experiment offers a pure and unique 
analysis of neutrino oscillations that is complementary 
to charged current-based appearance and disappearance 
searches. 
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